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Abstract 



We study the anomalous breaking of J7a(1) symmetry just above the QCD 
phase transition for zero and two flavors of quarks, using a staggered fermion, 
lattice discretization. The properties of the QCD phase transition are ex- 
pected to depend on the degree of Ua(1) symmetry breaking in the transition 
region. For the physical case of two flavors, we carry out extensive simula- 
tions on a 16 3 x 4 lattice, measuring a difference in susceptibilities which is 
sensitive to Ua{^-) symmetry and which avoids many of the staggered fermion 
discretization difficulties. The results suggest that anomalous effects are at 
or below the 15% level. 
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The breaking of classical Ua(1) chiral symmetry by quantum effects is a theoretical 
phenomenon of considerable physical importance which has direct implications on the order 
of the finite temperature, QCD phase transition. If we consider only the two light u and 
d quarks, anomalous symmetry breaking reduces the flavor symmetry from {7(2) x {7(2) to 
SU(2) x SU(2) x Ub{1) x Za{2). It is only this reduced symmetry that is consistent with 
second-order, critical behavior Although present lattice simulations suggest that the 
QCD phase transition is indeed second-order for two flavors, it is an important consistency 
check to establish the required anomalous £7a(1) symmetry breaking directly. In addition, 
the size of these symmetry breaking effects will determine the width in temperature and 
quark mass of the region showing universal critical behavior. 

The plasma phase of QCD is a particularly good place to study the axial anomaly. For 
T < T c , the dynamical breaking of chiral symmetry obscures the effects of the axial anomaly. 
For T > T c , chiral symmetry is restored and thermal Greens' functions are explicitly sym- 
metric under SU(Nf) x SU(Nf) transformations. We can then look directly for anomalous 
symmetry breaking by comparing Greens functions that are related by the anomalous {7a (1) 
symmetry. 

Such breaking of global {7a(1) symmetry, associated with a zero-momentum Ward iden- 
tity, is especially interesting, since at zero momentum the anomalous term in the chiral 
Ward identity becomes Nf times the topological charge u, a quantity which vanishes to all 
orders in conventional perturbation theory. The rj mass in QCD and the non-conservation 
of baryon number in the standard model are other examples in which such non-perturbative 
anomalous effects should occur 0. 

This sort of anomalous symmetry breaking can be understood from two perspectives: i) 
As the physical remnant of an ultraviolet ambiguity in the theory. Here modifications made 
to regulate the divergences present in the continuum gauge theory, necessarily introduce 
explicit chiral symmetry breaking whose effects remain visible, even on energy scales small 
compared to those of the regulator, ii) As resulting from an infrared singularity which 
permits chiral asymmetry to survive the chiral limit of vanishing fermion mass. While chiral 
symmetry breaking effects might naively be expected to be proportional to the explicit quark 
mass, m, the presence of infrared singular, 1/m behavior can allow such chirally asymmetric 
effects to remain even in the m — > limit. In a semiclassical calculation, such 1/m behavior 
appears for gauge backgrounds with fermion zero modes. These two views are related by 
the Atiyah-Singer index theorem |J. 

Most lattice calculations [Q explore the first approach, studying explicit chiral symmetry 
breaking inherent in the lattice regulation which should reduce to anomalous effects as the 
continuum limit is taken 0]. In this paper we explore the second approach, searching for 
anomalous asymmetries that arise from infrared singularities in the limit of small quark 
mass. For lattice calculations these two approaches are not equivalent since the Atiyah- 
Singer theorem applies only in the limit of an infinite number of degrees of freedom. 

It is important to pursue both methods. The first approach may over-estimate anomalous 
effects, confusing them with simple lattice artifacts which vanish in the continuum limit. The 
second approach may miss anomalous effects since infrared singularities are often softened 
by lattice effects, e.g. the zero-mode-shift of Smit and Vink ||, and become apparent only 
as a — > 0. 

The question of anomalous symmetry breaking above T c has now been studied by a 
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number of groups. For a general review see the article of Laermann 0. Preliminary versions 
of our results can be found in Ref . || , while an alternative calculation, also taking approach 



ii), can be found in Bernard, et al. ||. Finally, Kogut, et al. |T0| use a combination of 
methods examining signals for anomalous symmetry breaking of both types i) and ii) above. 

In this paper, we study both zero- and two-flavor QCD, just above T c . For Nf = 0, 
anomalous effects are expected in the chiral condensate, (qq), while for Nf = 2 we must 
examine a more infra-red singular, L^(l)-noninvariant quantity, here a difference of iso- 
vector susceptibilities which we refer to as to = xp — Xs where 

XP = 2^/ rf4xrf4 ^^7\(^)x^7 5 x(?/)) (!) 

for space-time volume Q and flavor generator tK The scalar susceptibility xs is defined 
similarly, by omitting the internal ij 5 factors. 

We adopt the staggered fermion, lattice discretization. This is the approach used most 
successfully to date in finite temperature, lattice QCD studies. While the quantity (qq) can 
be calculated directly using this formalism, the more physically interesting u = xp ~ Xs 
can not. Direct definitions of xp an d Xs using staggered fermions will necessarily intro- 
duce ambiguities, with potentially large lattice artifacts obscuring the anomalous effects of 
interest . 

Here, we take an indirect approach, expressing uj, in the continuum, as a spectral integral 
whose singular behavior as m — > gives rise to anomalous symmetry breaking. We then 
demonstrate that this spectral integral can be directly evaluated using staggered lattice 
fermions and use this result to provide a lattice calculation of u. 

Consider the spectral representations: 

M = -tmc r dx!%£™>.\ (2a) 

JO A 2 + ITL^ m ( =m 

U = 4m L d V + m')» - (2b) 
Here p(X,g 2 ,m) is the average density of Dirac eigenvalues A. The first formula is due 



to Banks and Casher JTI] and the second is derived in a similar fashion. In Eq. ^a] we 
distinguish the fermion mass that appears in the fermion line attached to q and q, m^, 
from that entering through the fermion determinant, m. The factors of m or in the 
numerators of Eq. |2| reflect the chiral symmetry breaking character of (qq) and uo. However, 
an anomalous, small-mass limit can result if the integral over A is sufficiently singular for 
small A. 

Now let us investigate what might be expected for these quantities in continuum QCD. 
For T > T c , the small mass limit of (qq) and to in the continuum theory can be analyzed for 
both the case of very small volume and in the limit of infinite volume. For finite volume, the 
Dirac spectrum will be discrete for each gauge configuration in the path integral. The only 
non-zero contributions to either (qq) or u as m — >• will come from gauge configurations 
with at least one exact Dirac zero mode. In very small volumes, these zero modes can be 
predicted semiclassically and give the anomalous, small-mass behaviors: (qq) ~ 1/m, for 
Nf = and u> ~ constant, for Nf = 2. 



3 



The case where V — > oo first is more interesting and can be analyzed using the methods 
of Leutwyler and Smilga [Q. Above T c , there are no massless modes so the free energy 
should be proportional to the volume and analytic in the fermion mass: 

Z « expfi{F + FatrM+M + Gre{e i9 detM}}, (3) 

where M is the complex fermion mass matrix and 9 the usual theta parameter. Defining 
the topological susceptibility as Xtop — 

one easily derives Xtop — 

ttGm N f for 

9 = and M = ml, a multiple of the identity. 

We can similarly obtain expressions for (qq) and u: 

(ft) = --^-^-lnZ(M) = -2F 2 m-Gm N f- 1 (4a) 
JSf \l dm 

oj = ±{-^ - -^}lnZ(M) = 2Gm N f-* (4b) 
i I ami om\ 



where in Eq. |4a| we have divided by Nf to define (qq) as coming from a single fermion specie 
while in Eq. |b| we have used a complex M = ml + {mP x + i m\ )r J '. 

If we make the possibly reasonable assumption that the quenched value of (qq) can be 
obtained as the Nf — > limit of Eq. [la], then we can combine Eq. |] with the formula for 
Xtop to obtain: 



Nf = : (gg) = ^ JL 

Nr = 2 : = 2 xtop _ congt 



The last relation is of particular interest, implying that above T c the quantity uj provides 
an alternative measure of the topological susceptibility. As is shown below, uj can be easily 
determined using lattice methods, without the normal difficulties of defining topological 
winding on a discrete lattice. 

We will now compare these continuum expectations with lattice calculations. Because 



of the remnant chiral symmetry of staggered fermions, Eq. |2a] is also valid on the lattice, 
allowing us to relate (xx) an d p, where x is the single component, staggered fermion field. 
Viewing (qq) as a function of m and m^, we can express uj as a function of (qq) and then 
use this continuum result to define u on the lattice: 

" = ~(XX) + ^(Xx)k- (6) 

where these two terms correspond precisely to the terms in the difference uj = xp — Xs- in 
the remainder of this paper we quote values of uj, xp and xs normalized according to Eq. |6| 
where (Xx) is normalized to behave as 1/m in the large mass limit. 

First consider Nf = 0. In Fig. [1| we show (Xx) f° r two distinct phases distinguished by 
the complex phase of the Wilson line, (W), computed at (3 = 5.71, just above (3 C = 5.6925. 
(Recall the Wilson line, W, is the volume average of the trace of the ordered product of 
link variables along a line in the time direction.) For the case where (W) is real, we see 
the power law ~ m a76 for both 16 3 and 32 3 volumes suggesting this power law description 
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holds in the infinite volume limit. We see no sign of the anomalous 1/m behavior in (xx) 
expected in the continuum. 

For the case of complex (W), we see an unexpected spontaneous breaking of chiral sym- 
metry above T c , with (xx) approaching a constant as m decreases. The eventual decrease 
in (xx) for very small m < m m ; n is the normal finite-volume behavior expected with spon- 
taneous symmetry breaking, with m m i n (xx)V/T ~ 1 for both volumes. 

Next we examine uj and the more physical case of two flavors, at f3 = 5.3, just above (3 C 
(recall j3 c ~ 5.265 for N t = 4 and ma=0.01), on a 16 3 x 4 lattice for five different values 
of the dynamical quark mass. The results are summarized in Table | and plotted in Fig|| 
This figure shows the chiral condensate, (xx) approaching zero linearly as is expected for 
(3 > /3 C . Likewise, xp shows the expected regular, constant behavior as m — > 0. However, 
rather than showing the anomalous behavior, uj ~ ujq + uJ2fn 2 , expected from Eq. [|, Fig.|2| 
suggests a nearly linear uj as m goes to zero. 

Four fitted curves are also shown in Fig. ^. The two linear fits to (xx) and uj have a 
X 2 /dof of 2.2 and 2.7 respectively. Both of these fits are constrained to vanish at m = 0. If 
that constraint is dropped for the uj fit, the intercept moves upward slightly to 0.15(5) and 
the x 2 /dof falls to 0.34. A fit to the expected form ujq + u^m 2 is worse, with a x 2 /dof of 
3.4. 

While these results are consistent with those reported by Bernard et al. ||, our conclu- 
sions are different. That calculation examines a smaller lattice spacing than considered here 
but with larger statistical errors. Their analysis adopts the quadratic small-mass dependence 
for uj. However, such quadratic behavior is only guaranteed theoretically in the unphysical 
limit where m vanishes at fixed lattice spacing and needs to be established numerically for 
the case of interest. 

From Fig. one observes that for quark masses in the range 0.01 < m < 0.025, our 
results are consistent with an unusual but non-anomalous, linear behavior uj ~ m. At our 
smallest mass, 0.005, uj is significantly higher than such a linear extrapolation, suggesting 
that anomalous effects may be emerging. However, such effects are clearly quite small and 
occur for quark masses that are below those used in present studies of QCD thermodynamics, 
suggesting little connection between this anomalous behavior and the observed second-order 
QCD phase transition. 

In order to describe the physical size of a possible non-zero value of uj\ m= o, we must 
address the potential cut-off dependence of the quantities being discussed. While a thorough 
analysis of this question lies beyond the scope of the present paper |T3|] , there are two issues 



that are important to recognize. First, the m = intercept of uj requires the same ln(a)- 
dependent, multiplicative renormalization as the inverse square of the quark mass m, as is 
suggested by Eq. j|. We will ignore such a factor for our present rough estimate, since this 
factor should be of order 1 for current lattice spacings. 

Since uj is the difference of xp and xs, h is natural to compare uj to either of these 
quantities. However, both quantities contain a 1/a 2 piece when evaluated in physical units. 
Thus, we choose to compare wtoa quantity we will call uj, obtained as the difference between 
Xp evaluated at (3 = 5.3 and xs evaluated at (3 = 5.245, just below the transition: 



~ _ bod i _ d (xx) I m 
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Such a subtraction removes the unwanted, quadratically divergent 1/a 2 term at tree level 
and leaves an expression finite up to an 0(1), ln(a)-dependent multiplicative factor and a 
much smaller 1/a 2 term suppressed by the factor (5.3 — 5.245). We find uo/uj ~ 15%. 

Within the expected critical region, the light modes (7?, a) should be much less massive 
than the non-universal degrees of freedom suggesting uj/uj ~ 1, not the ~ 0.15 observed 
here. Thus, our results suggest that 0(4) critical behavior should not be seen in Nt — 4 
thermodynamics at least for \/3 — (3 C \ ~ 0.03. 

In conclusion, we have numerically studied anomalous symmetry breaking by examin- 
ing quantities whose anomalous behavior comes directly from infrared effects. Given the 
relatively coarse lattice spacing in our simulations a ~1 Fermi, our failure to find such ef- 
fects above the 15% level is far from conclusive evidence that such effects are suppressed 
in Nature |IJJ. However, this represents a first step in a systematic lattice calculation of 



such phenomena and must be followed by more demanding calculations on finer lattices and 
calculations using fermion formulations with improved chiral properties. 

We thank Edward Shuryak and Andre Smilga for helpful discussions and Yubing Luo for 
assistance. 
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TABLES 



TABLE I. Our (3 = 5.3, 16 3 x 4 results for two flavors of dynamical quark with mass m. Run 
length is the number of time units in the hybrid, 'R'-algorithm evolution after 200 time units were 
discarded. These quantities are normalized in a manner consistent with Eq. |6| with \X defined so 
that it behaves as 1/ma for large ma. 



ma run length {xx) Xs u 

0.005 4464 0.02256(21) 3.945(43) 0.559(37) 

0.01 2550 0.04374(50) 3.369(68) 0.932(56) 

0.015 2600 0.06517(82) 3.010(52) 1.322(62) 

0.02 2992 0.0896(10) 2.697(34) 1.722(62) 

0.025 3072 0.1141(38) 2.38(11) 2.24(12) 
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FIG. 1. The chiral condensate (xx) plotted as a function of quark mass for a pure gauge 
calculation on 16 3 x 4 and 32 3 x 4 lattices. The real phase (closed points) is the most physical 
(det(D — m) is largest for this phase). No evidence is seen for the expected anomalous behavior, 
(xx) ~ m~ l as m — > 0. 
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FIG. 2. The quantity u, which directly measures anomalous symmetry breaking, plotted versus 
fermion mass, ma. Also shown are the chiral condensate (xx) an d the pseudoscalar susceptibility 
Xp- We studied a 16 3 x 4 lattice at /3 = 5.3, just above j3 c - 
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